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Abstract 



• Some basic difficulties and ambiguities on the way from the mea- 

sured momentum distributions to the inferred properties of the inter- 
^ | action regions are reviewed and discussed. 

^ 1 Introduction 

The Bose-Einstein correlations in multiple particle production processes are 
studied mainly in order to get information about the interaction regions, i.e. 
about the regions where the final state hadrons are created, and about the 
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hadron's last scattering surfaces which by definition are the boundaries in 
space-time of the interaction regions. There is a number of things one would 
like to know about an interaction region: its size and shape, its orientation 
with respect to the event axis and in the case of the heavy ion collisions 
also with respect to the impact parameter vector; the flows of matter, if 
any, within the region; the nature of the matter filling the interaction re- 
gion, its equation of state and its phase transitions, if any. This list could 
be made longer. One could try to obtain this information either from first 
principles, or using phenomenological models, or in some model independent 
way from the data. The approach from first principles is for the moment 
too difficult. In the present paper we show some of the difficulties in deriv- 
ing information about the interaction region directly from the data without 
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model assumptions. The conclusion is that little reliable information can be 
derived without a trustworthy model. The data in practice means momen- 
tum distributions: single particle, two-particle, three-particle etc. There is 
some additional information, e.g. every experiment gives some information 
about the positions of the interaction vertices, but when we are interested 
in distances of the order of tens of fermis or less this additional information 
does not contribute significantly. The fc-particle momentum distribution is 
given by the diagonal elements of the fc-particle density matrix in the mo- 
mentum representation. Most models assume that this can be expressed as 
a symmetrized product of single particle density matrices pQ: 

k 

p(pi, , p k , Pi, , p fc ) = II Pi(Pi5 Ppj)' ( x ) 

P j=l 

where the summation is over all the permutations of the momenta p 1; . . . , p k . 
Thus, e.g. the single particle momentum distribution is pi(p; p) and for the 
two-particle distribution we get the well-known formula 

p(pi,p 2 ;pi>P2) = pi(pi;pi)pi(p 2 ;p2) + bi(Pi; p 2 )l 2 - ( 2 ) 
2 Emission function 

A popular strategy when studying the Bose-Einstein correlations is to evalu- 
ate the emission function S(K, X) (cf. 2\ and references given there) related 
to the single particle density matrix by the formula 

p 1 (p 1 ;p 2 )= JdXS(K,Xyi x . (3) 
Here K, q and X are four-vectors defined by 

K = ^{Pi+P2)] q = Pi-P2] X = ^(x 1 + x 2 ). (4) 

Note that Kq = 0. It is usual to interpret S(X, K), for given K and x = t, 
in the spirit of the interpretation of the Wigner function, as the phase space 
distribution of particles. Then S(K, X) gives directly all the information 
about the geometry and the evolution in time of the interaction region. 
Moreover, it gives the momentum distribution of the produced particles and 
information about the momentum-position correlations. The emission func- 
tion can be calculated when the production amplitudes are known [3], jl]. 
Let us consider now the prospects for calculating the emission function from 
the data using formula (jHJ). 



2 



On both sides of equation (jSJ) the four-vector K has the same, fixed value. 
The right-hand-side is a Fourier transform of S changing the variables into 
the variables q^. This Fourier transformation, however, cannot be inverted, 
because the left hand side is known only for the four-vectors q satisfying 
the condition Kq = 0, i.e. for only one value of go at each q. This implies 
that for any given density matrix p\ there is an infinity of different emission 
functions satisfying relation (jHJ). They correspond to the infinite variety of 
possible continuations of function pipxWi) t° unphysical values of go- The 
data without additional information from theory, or from a model, are unable 
to tell us which of these emission functions is the good one. In order to see 
how this affects the conclusions about the interaction region let us consider 
the following example. 

We replace equation (j3J) by the equation 



where c is an arbitrary real constant. This is, of course, legitimate, because 
besides renaming the integration variable X' we have just added to the ex- 
ponent icKq = 0. Changing the integration variable to X = X' — cK we 

get 



Thus S(K, X + cK) is another solution of equation Q for the same density 
matrix p\. The space-time distribution of the sources (more precisely of the 
particles produced by the sources) is obtained by integrating the emission 
function over K. Suppose now that for c = there are no X — K corre- 
lations, i.e. for every K the space time region occupied by the interaction 
region is the same. Then the K integration does not affect the space time 
distribution — one can say that the interaction regions for all the K-s are 
piled on top of each other. When c becomes different from zero the interac- 
tion regions corresponding to the different K values get different shifts. As 
a result the size of the overall space time region occupied by the interaction 
region increases. With sufficiently large c it can be made as large as one 
wishes. For instance, one could make it the size of a football, or living for 
an hour. Such extreme choices of c can be eliminated using what we know 
about the vertex positions in space-time, but using momentum distributions 
only, they are just as good as c = 0! 

The obvious question is: what information besides the momentum mea- 
surements is necessary to obtain reasonably accurate descriptions of the in- 
teraction regions? One could try to evade the problem by saying: just forget 
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about the emission functions and your problem will disappear. We will show 
in the next section that this is not a satisfactory answer. 

3 Density matrix and Wigner function 

Let us assume that all the hadrons have been produced in their final states 
(scattered for the last time) instantaneously and simultaneously at some time 
which we may chose as t = 0. Then, for any t > 0, all these hadrons can be 
described by standard, time independent (in the interaction representation) 
density matrices. According to assumption all these density matrices 
can be expressed in terms of the single particle density matrix p\. Our key 
observation [5 j is that replacing p\ by 



where a(p) is any real- valued function of momentum, does not change the 
momentum distributions. This is obvious for k — 1 and k = 2. In general, 
however, on the right-hand side of equation (0), in every monomial every 
momentum p i occurs exactly once as the first argument of p\ and exactly 
once as the second argument of p 1 . Thus, when pi gets replaced by pi a , it 
brings two factors e^P^ and e -10 ^ 1 ^ which cancel. 

As a corollary let us note that this makes the definition of the emission 
function, if obtained from the data using equation (jSJ), even more ambigu- 
ous. Not only for given pi there is an infinity of different solutions for S, but 
moreover p 1 can be replaced by any of the infinity of functions p\ a without 
affecting the fits to the experimental data. Using the density matrix formal- 
ism instead of the emission functions, we have the second ambiguity, but not 
the first. Let us see how it affects the inferences concerning the interaction 
regions. 

Let us make the simplifying assumption that 



Using the general formulae for the diagonal elements of the density matrix 
in the coordinate representation and for the Wigner function in terms of the 
elements of the density matrix in the momentum representation: 




(7) 






(10) 



(9) 
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one gets 



pi(x;x) 
W(K,X) 



1 



1 



x 



(2vri?A) 3 



cxp 



exp 



2R 2 
K 2 



X 2 



2A 2 2R 2 



(11) 
(12) 



The diagonal elements of the density matrix in the coordinate representation 
yield the distribution of particles in space at t — (keep in mind that we 
are using the interaction representation). This is a spherically symmetric 
Gaussian with the root mean square width in every coordinate equal R. 
The Wigner function is usually interpreted as an approximation to the phase 
space distribution. In our case it confirms the result for the space distribution 
and, moreover, gives the information that the momentum distribution is 
uncorrelated to the space distribution and is also Gaussian with the root 
mean square width A. 

Let us now replace the Gaussian p 1 by the corresponding p la with 



-cp 2 , 



(13) 



where c is an arbitrary real number. A simple calculation gives 



^2ir(R 2 + c 2 A 2 )' 



Pi Q (x; x) = 

v 2 



(27T.RA) 3 



cxp 



• exp 



K 

"2 A 2 



2(i? 2 + c 2 A 2 ) 

W a (K,X) = 
(X + c K) 
2R 2 



(14) 



(15) 



The first formula shows that the distribution in space remains spherically 
symmetric and Gaussian, but the mean square radius can be any number 
not less than R. The second equation shows that momentum-position corre- 
lations have been introduced and that the smearing of the interaction region 
is due to a mechanism similar to that described in the preceding section: the 
interaction regions for particles with different values of the momentum K 
shift with respect to each other. 

Note that, if we have a model which predicts unambiguously the den- 
sity matrix pi, there is no problem. The problem arises if one insists that 
the conclusions about the interaction region should be deduced from the 
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data without additional assumptions. An interesting question, however, is: 
suppose that we have two models based on completely different physical as- 
sumptions, one giving a matrix p\ which agrees with experiment, and the 
other a corresponding matrix p\ a which, as follows from the preceding dis- 
cussion, agrees with experiment just as well; how can we tell which of the 
two models is the good one? Experimental data on momentum distributions 
cannot help however precise they are. The standard advice to look at more 
particle distributions is also useless. 

4 Distribution of distance between points where 
identical particles are produced 

Let us consider a pair of identical particles with total momentum 2K and 
denote the distance between their production points by x. An important 
strategy is to study the distribution of the distances (x), i.e. 5j^(x) [2], 
0. The relation of the function Sj^ to the data is 



dN/d Pl dp 2 f 3 r 2 -I 

dN/d Pl dN/d P2 = r XiS k(x ) [i<M x )i - x \ > ( 16 ) 

where the primed quantities are defined in the center of mass frame of the 
pair and 0q>(x') is the two-particle wave function of the pair in its center 
of mass frame. When two particles are far from each other their relative 
momentum is well defined and equals q'. In particular for free particles the 
content of the square bracket equals cos(q' • x'), but one can include final 
state interactions by using more complicated wave functions. 

The solution of equation (jlfij) for the free particle case can be easily done 
by inverting the Fourier transformation. For more complicated cases one has 
to use more advanced methods known as imaging [7j, jH], [H], but there is 
no difficulty of principle. In this approach function Sj^ is unambiguously 
given by the data. This is a nice change after the two methods described 
in the preceding two sections and, therefore, the method is rapidly gaining 
popularity. 

Let us note, however, that the left hand side of equality (fTBj) . in approx- 
imation (Q) is 

dN/dp x dp 2 _ \pi(pi,P2)\ 2 +1 ^ 



dN/dp 1 dN/dp 2 pi(pi, Pi)pi(p 2 ; P2) 

The right hand side is clearly invariant with respect to the transformation 
discussed in Section 3. Therefore, there is an ambiguity concerning the con- 
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elusions about the interaction region as discussed there. We infer that func- 
tion can be evaluated unambiguously from the data, but its implications 
for the interaction region are highly ambiguous. Thus, again there is no sim- 
ple way from the data to the information about the interaction region. It 
is plausible that each of the three approaches discussed here, and certainly 
many others, share essentially the same difficulty, which just changes its form 
when the formalism is being changed. 



5 Conclusions 

The analysis presented in the previous three sections strongly suggests that 
little information about the interaction region can be obtained from the 
data on the momentum distributions without making additional assump- 
tions. Usually this problem is solved by using a model. Every model contains 
assumptions impossible to prove by just using the data on the momentum 
distributions. Within the model these assumptions could be justified by ar- 
guments of beauty, simplicity, or by invoking some general principles. This 
rises the following important problem, however. 

Suppose that two models differ only in these additional assumption and 
thus they give exactly the same predictions for all the momentum distribu- 
tions. Suppose further that they correspond to widely different and con- 
flicting pictures of the interaction region and of what happens there. How 
should one decide which model is the realistic one? This is not an academic 
problem. It has been already noticed that widely different model can fit the 
same data, cf. e.g. [10], [TTj . 
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